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Abstract 

The high strain rate mechanical response of thin polymer plates 
has been studied using a modified shock tube. Diagnostics include 
the pressure-time history of the incident and reflected pulses and the 
use of digital image correlation (DIC) techniques to extract the time- 
history of the out-of-plane displacement distribution. Additionally, 
finite element models have been developed to understand the plate 
response and to validate and modify plate material constitutive mod- 
els that have been proposed. 

Key Words: Shock Loading , High Strain rate , Shock Impact , Material 
Characterization 


1 Introduction 

The study of the response of materials subject to blast loading is impera- 
tive to the understanding of how to build blast resistant materials. Until 
recently the only way to test the materials was through actual blast tests 
on large size panels, that are costly. However recent developments have 
used a modification to a shock tube to simulate blast conditions, provid- 
ing an alternative way to investigate material behavior subjected to blast 
loading [1]. 

The shock tube is a well understood instrument, with much of the the- 
ory and validation available in the early literature [2]. The corresponding 
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fluid dynamic equations can be found in [3]. In the present work, a mod- 
ified version of the shock tube is used to understand how a shock wave 
deforms an initially flat panel. By removing the very thick end cap of 
the driven tube and replacing it with a relatively thin panel specimen, it 
becomes possible to subject the panel to a very large pressure pulse act- 
ing over a very short-time duration. This method of loading is a very 
recent development and the literature reflects the novelty of the test at 
this point [1]. Many of the initial tests conducted using a shock tube were 
only minimally instrumented and subsequently post-shock tests were con- 
ducted to determine the compressive strength or other material property 
degradation [4,5]. Additional shock tube tests have been performed on 
composite plates [6-9], however a lack of data during the deflection and 
straining of the plate makes it difficult to evaluate the properties and re- 
sponse. 

Shock tube studies in conjunction with the development of finite el- 
ement models have been previously investigated. Initial studies focused 
on validating an anisotropic damage evolution model for metal plates [10]. 
Another study attempted to match experimental data of viscoplastic plate 
vibrations caused by the shock wave impact with finite element simula- 
tions [11]. Ultimately the goal would be to blend both problems, i.e that 
of the shock interaction with the plate response, which becomes a fluid- 
structure interaction problem. Initial attempts to combine both a fluid dy- 
namics code with a finite element code have been undertaken, with very 
primitive results [12]. On the other hand if the timescale associated with 
the shock wave loading is significantly smaller than the plate response 
time, then the shock loading approximates to a steep pressure step pulse 
and the ensuing response, if accurately measured, can be used to validate 
and refine the constitutive models for a variety of rate-dependent materi- 
als. 

In this paper, results from shock wave loading of epoxy panels are pre- 
sented. The time-resolved panel response, which is accurately measured 
using digital image correlation (DIC) techniques, is studied in conjunction 
with a finite element simulation of the response. This provides a means to 
validate a rate-dependent constitutive model for the epoxy material. 
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2 Method 


Figure 1, shows a schematic of the shock tubes in the Composite Struc- 
tures Laboratory (CSL) at the University of Michigan. The driver section 
is loaded with a high pressure gas which is separated from the driven 
section by a diaphragm designed to burst at a desired driver pressure. 
When the diaphragm ruptures, a shock wave is formed which travels 
down the tube at a much faster rate than the pressure can equilibrate. 
When the shock wave hits the panel specimen, the specimen surface expe- 
riences a sudden rise in the pressure and temperature. The shock tube is 
equipped with pressure transducers (see figure 2), along the driven tube. 
The pressure transducer mounted on the tube is located 9 in.(0.23 m) from 
the location of the specimen. The pressure transducers used were made 
by PCBPiezotronics (model number 113A22) which have a measurement 
range between 0-5000psi (0-34.5 MPa ) and a resonant frequency greater 
than 500kHz. They are connected to a 4 channel signal conditioner (model 
number 482 A22). The signal from the signal conditioner then enters the 
oscilloscope. The driven section always remains at atmospheric pressure. 

A complex viewing port system for providing lighting, using flash- 
bulbs, and access for high speed cameras to obtain images to perform the 
DIC measurements is also present. A schematic of the shock tube hard- 
ware incorporating the camera, that was used to obtain the DIC measure- 
ments, is shown in figure 3. Prior to testing, the back surface of the speci- 
men is sprayed with a speckle pattern that is suitable for subsequent DIC 
measurements. 

The DIC system provides full-field displacement measurements on the 
back surface of the panel at pre-selected instances as the deformation pro- 
gresses. The images are captured using two synchronized Photron SA.l 
cameras, that capture images at a rate of one hundred thousand frames per 
second at a resolution of 320 x 128 pixels. The displacement information 
can be used to determine the strain components; e x ,fY, 7 xy as a function 
of position and time on the back surface, and also the time rate of change 
of these components at different locations on the panel back surface. This 
is very useful since different areas of the panel experience different levels 
of strain and rates of straining, depending on the spatial location. 

The accuracy of the measurements can be related to the facet size in 
the DIC software. The resolution of the measurement is 0.16m.( 4 mm). 
Although the step size is 0.0094m. (0.238mm) the size of the computation 
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Figure 1: Schematic of the shock tube setup. 
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Figure 2: Pressure transducer locations on tube and the plate (specimen). 
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Figure 3: DIC setup for Shock tube Measurements. 
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box will determine the spatial resolution. The random speckle pattern uses 
a dot size that is roughly 5-7 pixels in diameter, as recommended by the 
software manufacturer. The software used in this analysis was ARAMIS, 
provided by GOM. 

The epoxy material was provided to the University of Michigan by 
NASA Glenn research center. No details on the type of material nor the 
expected mechanical properties, etc. were provided. It was simply pre- 
sented as a matrix that would be used in a composite. 

Additional tests were performed with a steel plate that was 0.75 in.(19 
mm) thick, with mounting holes situated as specified by the manufacturer, 
at the center, half the radius and edge, to determine the pressure experi- 
enced on the surface of the plate. Test results showed that the shock pres- 
sure on the plate was uniform, allowing us to assume that the shock wave 
was normal to the surface. 


3 Experimental Results 

A series of tests were conducted using different thickness epoxy plates. 
For each thickness test, identical shock tube diaphragms were used so that 
the burst pressure, and hence the pressure pulse on the specimen, would 
be nearly the same for materials of the same thickness. The time history of 
the maximum out-of-plane displacement for a representative test can be 
seen in figure 4 along with the diaphragm burst pressure. Also indicated 
is the maximum pressure imposed by the specimen. These results were 
obtained using the DIC software, and the results of a single frame can be 
seen in figure 5, which corresponds to a known pressure history of the 
thickest panel. 

In addition to the panel tests performed, another test was run with a 
"rigid" plate mounted with pressure transducers to record the pressure 
that was imposed on the impacted plate surface. This pressure-time his- 
tory is assumed to correspond to those experienced by the epoxy panel 
specimen. This assumption is valid if the epoxy panel response time is sig- 
nificantly larger than the rise time of the shockwave pressure pulse. The 
two tests ('rigid' panel and epoxy panel) were performed to experimen- 
tally determine the effect of panel deflection on the pressure felt by the 
panel (i.e. the fluid-structure interaction issue). This information is very 
important in the subsequent finite element modeling, since one needs to 
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Figure 4: Comparison of maximum out of plane displacement for various 
specimen thicknesses. The pressure denotes the corresponding diaphragm 
burst pressure for the test, and the maximum pressure on the specimen. 



[in] 


- 0.00643 

- 0.01500 

- 0.02250 

- 0.03000 

- 0.03750 


- - 0.04500 
-- 0. 05250 
(- 0.06101 


Figure 5: Contour plot of out of plane displacement deformation for a pure 
epoxy panel. 
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know the input pressure time-history. Figure 6(a) shows the two pressure- 
time histories (that have been shifted in time) corresponding to the rigid 
plate and the thinner epoxy plate as measured. The pressure time history 
measured on the tube is shown in figure 6(b) along with the pressure mea- 
sured on the surface of the plate. The results show that the pressure-time 
histories are very similar, even on the deformable panel. One clear reason 
for this is the time scales involved in the problem. The pressure pulse in- 
crease occurs over a time duration of 0.5 /is, while the panel deflection, in 
response to this pressure pulse, occurs over a time duration of 250 /is, as 
can be seen in figure 6(c). In summary, the reflected pressure that is mea- 
sured by the tube pressure transducer can be used as the pressure pulse 
experienced by the panel specimen. It should be noted that this does not 
account for fluid-structure interaction that can occur from a rapidly de- 
forming material or interaction with a specimen that fails. 

A comparison of the pressure-time history, to the center deflection-time 
history is shown in figure 6(c). Since the center deflection- time history 
shows a clear separation in scale, the center plate pressure now appears 
very similar to a square wave. This result is important in simulations be- 
cause the shock pulse can be modeled as a square wave to determine what 
pressure level is needed to induce the desired strains and stresses in a 
given plate specimen. This calculated pressure can then be related back to 
a needed burst pressure to produce the required pressure pulse. 

With an understanding of the pressure-time histories, the pressure im- 
posed on each panel during testing can be evaluated. The results are 
shown in table 1 with the corresponding reflected shock pressure being 
measured from the pressure transducer mounted on the shock tube. The 
variation in specimen pressure has to do with how the shock wave is re- 
flected off the specimen. If the shock has enough time to compress the 
specimen the pressure will build however, if the specimen fails or has 
too great of a deflection, the shock wave will not be completely reflected. 
Since the thinner specimens failed with larger deflections the strength of 
the shock wave will not be as strong as if a rigid plate had reflected the 
shock wave. 
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(b) Pressure Transducer in tube vs Pressure 
transducer on surface of rigid plate 



(c) Pressure on rigid plate and center deflec- 
tion of epoxy panel 

Figure 6: Pressure- time histories and center deflection (Epoxy Panel) of 
shock tube test. Time has been shifted to show relative data comparison. 
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Table 1: Burst and shock pressures for various thicknesses. 


Thickness 

(in.) 

Diaphragm 

Burst 

Pressure 

(psi) 

Shock 

Pressure 

(psi) 

0.060 

1230 

140 

0.125 

630 

120 

0.192 

1215 

225 

0.242 

990 

240 


4 Computational Model 

A finite element (FE) model was developed to understand the pulse-loaded 
panel response. Initially, a homogeneous isotropic plate model using non- 
linear geometry was created in ABAQUS version 6.8, to simulate the epoxy 
panel response, using an explicit solver. The measured pressure-time his- 
tory is used as a "square wave" input loading. This assumption does not 
allow for fluid-structure interaction and the propagation of further shock 
waves. The epoxy material was modeled using J2 small strain theory of 
plasticity with isotropic hardening and a von Mises yield criterion, where 
the material properties are to be determined through inverse modeling, by 
correlation with the test data. A circular 3 inch specimen corresponding 
to the specimen dimensions used in the shock tests, was analyzed using 
the FE method. The panel mesh contains 458 elements with 254 nodes, 
with S3R shell elements being used. The mesh size was determined using 
a convergence study, the results of which are shown in figure 7 

For materials which have a rate-independent modulus the Cowper- 
Symonds overstress power law is used to approximate the effects of strain 
rate on yield stress. 


e pl = D(R - l) n , (1) 

where D and n are material constants and R is the yield stress ratio, 
Since the rate dependence of the epoxy is an unknown, an initial attempt 
was made using the rate-dependent properties of a similar polymer ma- 
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Figure 7: Finite Element Convergence Study 


terial. Split Hopkinson pressure bar test results for the rate dependence 
of SC-15, a thermoset polymer matrix material, were used as a guide for 
an "initial guess" of the rate-dependent properties. SC-15 epoxy is a low- 
viscosity, two-phased toughened epoxy resin system consisting of phase A 
(resin mixture of diglycidylether epoxy toughener) and phase B (hardener 
mixture of cycloaliphaic amine poluoxylalkylamine) [13]. The assumed 
rate dependence for SC-15 is plotted as a solid line in figure 8, where the 
circles represent the experimental data. 

Due to the complexity of the constitutive model many parameters need 
to be determined; these are, modulus, yield stress ratio, power law expo- 
nent, and material damping. The known inputs can be seen in table 2. In a 
brittle epoxy sample, a test that subjects the panel to failure provides very 
little information about the yield and plasticity parameters. From such a 
test, about all that can be determined are the initial modulus and failure 
strength. Therefore, if a panel test is conducted at a suitable burst pres- 
sure below the plate to fracture level, the resulting oscillatory out-of-plane 
panel response provides information that can be used to determine several 
of the material parameters. 

Through the use of inverse modeling the following parameters for the 
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Figure 8: e vs. ^ for SC-15 Matrix 
Table 2: Sample Test Known Parameters 


Thickness 

0.242 in. 

Density 

1.14e — 4 

in a 


epoxy material have been determined (see table 3), by matching the time 
history of events. With the use of the parameters shown in table 3, a com- 
parison between experiment and simulation for the panel out-of-plane de- 
flection can be made and is shown in figure 9. A consistent error between 
the computed and measured displacement amplitude history is seen, even 
though the frequency of response is matched almost exactly. Reviewing 
the hardware configuration of the DIC system, it becomes apparent that 
an inherent error was inadvertently introduced into the system. Figure 
10 shows that the optical path can shift, due to any shift in the position 
of the lexan panel window, causing a parallax error. When the system 
was initially calibrated the containment chamber was a distance 'd' away 
from the mount. However when the tests were run 'd' became zero as the 
system was clamped together. The cameras remained at the same fixed 
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Figure 9: Finite Element Results vs. DIC results 


distance from the point where the specimen is mounted. This shift will 
cause a change in the DIC systems measurements since the optical path 
shift introduces a systematic error into the initial image calibration. When 
the DIC system is initially calibrated, the distance 'd' must be held at the 
same value to be used later during the actual tests. Additional error could 
have been introduced if the assumed material properties from the SC-15 
data were inaccurate, since only one run was used to curve fit many dif- 
ferent properties. 


Table 3: Sample Test Parameters 


Modulus 

350, 000 

lbs 

in 2 

Gy 

3,250 

lbs 
in 2 

D 

4880 


n 

2.883 



In addition to displacements, the strain histories are also compared. 
Since strains are based on displacement data there will again be a mis- 
match due to the parallax error in terms of magnitude however, the gen- 
eral trend and contours should be the same. The results compare a section 
line through the center in the DIC results, which are seen in figures 11 
and 12. There is an excellent agreement between the FE predicted result 
and the DIC extracted result. The values of e xx show nearly identical con- 
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Figure 10: Difference in optical paths from movement of lexan panel. 


tour maps which is due to very little out-of-plane displacement 0.06m. 
occurred relative to the overall size of the specimen 3 in. therefore, the cal- 
culations for strains becomes more reliant on the 2D motion, which can 
be captured accurately with 1 camera. The DIC results show that there is 
some localized strain effects, which could be caused by imperfections in 
the material. 

With a focus on strain rate dependent material properties, a represen- 
tative point in the center of the panel is examined. The center of the panel 
produces a biaxial stress state, which implies that both e xx and e yy are the 
same value for an isotropic material. The strain-time history for the center 
point can be seen in figure 13. As can be seen there is very good agreement 
between the two results. From the strain- time histories, the strain rates are 
computed and shown in figure 14. As can be seen in figure 13, there is a 
considerable amount of noise in the DIC data. The data was filtered, using 
a ninth order polynomial curve fit, to eliminate much of the noise so that 
one can observe a general trend. There is an initial peak followed by a spi- 
ral pattern as the plate oscillates. To better see this, the strain components 
at the center, as predicted from the finite element simulation, are shown 
in figure 15. This plot shows that we have attained a medium strain rate 
(50-200 1, see figure 14) for the material. We have not achieved an ultra- 
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Figure 11: Comparison of e xx results at max displacement. 
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Figure 12: Comparison of e yy results at max displacement. 
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Figure 13: Strain time history of center point comparing FE vs. DIC 



Figure 14: Strain rate vs. strain history at the crown of the panel; compar- 
ison between the FE model and DIC results, which have been filtered. 

high rate (2000-3000 - s , see figure 8) as seen in the Hopkinson bar. The 
material will have ranges of constant strain rate and also oscillate about 
this strain rate. Figure 14 shows how the material will respond to a strain 
rate of 100^ staying relatively constant until it reaches a maximum strain, 
reversing and going to a strain rate of —75^, and then as the oscillations 
continue it will return to 75- and shift to -75- as the deflection decreases. 

S S 

From figure 15 it is observed that the shear strain is negligible, since its 
values and contributions are less than 1% of e x and e y . 

The remaining test results (where the panel fractured) were used to de- 
termine the stress level when failure occurred. Since this material is very 
brittle at high strain rates, there is minimal yielding. The test results were 


N AS A/TM— 20 10-216941 


17 


200 



Figure 15: Strain rate vs. strain history at the crown of the panel from FE 
model. 

compared to properties determined by the FE simulation model and used 
to back out the failure stress, based on the time that failure was recorded 
with the Photron cameras. It should be noted that the model is based 
on SHPB data from a separate material, along with the fact that the data 
was collected at much higher strain rates. Additionally, the different thick- 
ness of materials may have different material properties, based upon slight 
processing variations. Therefore, the model should only be used to note 
trends. 

A summary of the computed failure stresses, and the out-of-plane de- 
flections, as a function of panel thickness, is provided in table 4. Here the 
failure stress is produced from the FE model and the displacement and 
failure time are obtained from the DIC software. It is seen that at the very 
high rates, panel failure occurs in the brittle-elastic regime, as was evident 
in the failure mode of the material and the small amount of plasticity mea- 
sured in the model. 


5 Conclusion 

A methodology has been developed to determine the mechanical response 
of materials subjected to pressure pulse loading. Using high speed cam- 
eras, the out-of-plane deflection time history has been measured using the 
DIC method. This result, in conjunction with a finite element model and 
Hopkinson bar data can be used to inverse model the panel response and 
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Table 4: Failure Stress 


Thickness (in.) 

Failure Stress (psi) 

Time (/is) 

Z-Disp (in.) 

0.060 

7,635 

120 

0.152 

0.125 

6,860 

120 

0.131 

0.192 

6,599 

130 

0.0874 

0.242 

5,421 

130 

0.0587 


validate the rate-dependent properties of the panel. The pressure-time 
history is used as an input to the finite element simulation. It has been 
shown that the pressure transducer closest to the panel on the driven tube 
provides information regarding the loading of the specimen. In the shock 
tube tests, it is beneficial to stop the testing prior to failure, rather than hav- 
ing an oscillatory panel response, since that allows the determination of a 
larger number of material parameters. Tests to failure provide additional 
information on rate-dependent material strength and mode of failure. 
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